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i 1 Abstract 

The phase diagram of a fluid confined between a planar and a conical walls modelling 
the atomic force microscope geometry displays transition between two phases, one with a 
liquid bridge connecting the two walls of the microscope, and the other without bridge. The 
structure of the corresponding coexistence line is determined and its dependence on the value 
of the line tension coefficient is discussed. 

a ; 

1 Introduction 

The behavior of fluids confined in nanosized structures depends sensitively on the properties of 
enclosing walls [TJ [2] . The equilibrium paradigm of such influence is provided by the capillary 
condensation phenomenon in a slit. Despite thermodynamic conditions favouring the gas 
phase in the bulk, the liquid phase may fill the space between the slit walls [3]-[7]. The shift 
of the gas- liquid coexistence line /i = fJ>h(T) in the slit with respect to its position in the bulk 
■ system \x = Ho(T) is described by the Kelvin law: ^o(T) — fJ,h(T) ocl/h, where h denotes the 

distance between the walls [8]. 

When the walls confining the fluid are non-planar and are like in the AFM geometry then 
one may expect formation of liquid bridge linking the opposite walls 9 I15| at thermodynamic 
conditions that favor gas phase in the bulk. The size and the shape of the liquid bridge depend 
on system's geometry as well as on the thermodynamic state specified, say, by temperature 
and chemical potential. There exists an additional attractive capillary force acting between 
O ■ the walls linked by the bridge which must be taken into account during the Atomic Force 

, Microscope measurements [16]- [20]. Formation of the liquid bridge is also exploited in dip-pen 

nanolithography [21]- [24]. There the existence of the bridge enables the flow of particular type 
of molecules from the tip of the AFM onto the planar substrate. The size of produced patterns 
depends on the geometry of the bridge which itself is controlled by the thermodynamic state 
of the system and its walls' geometry. These two examples already show that the knowledge of 
precise thermodynamic and geometrical conditions for the existence of liquid bridges together 
with their morphological properties is of crucial importance. 

In this paper we investigate the morphological phase transitions taking place in 3d fluid 
confined by the walls resembling the AFM geometry. They consist of a planar substrate and 
the surface of a cone modelling the microscope's tip. Our study is based on a macroscopic 
approach in the grand canonical ensemble with the purpose of constructing the relevant phase 
diagram and determining the possible shapes of the liquid-like bridge spanned between the 
planar substrate and the conical tip. In particular, we discuss the influence of the line tension 
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on the structure of the coexistence curve and examine the special cases related to the filling 
transition. 

2 Shape of meniscus 

We consider fluid confined between two inert substrates, see FiglH The surface of the lower 
substrate (1) is an infinite plane z = and the upper substrate (2) forms an infinite cone 
characterized by the opening angle it — 2ip. The distance between the cone's apex and the 
plane z = is denoted by h; the parameters h and (p completely characterize the system's 
geometry which has cylindrical symmetry. In cylindrical coordinates the cone is described by 
equation r = a(z) = (z — h) cot <p, where r = \J x 2 + y 2 and z > h. The system under study 
resembles the AFM-like geometry in which the conical substrate (2) plays the role of micro- 
scope's tip. In the limit tp = 0, h ^ the geometry becomes that of a slit with two parallel 
walls separated by the distance h. In another limiting case corresponding to ip ^ and h = 
the confining walls can be obtained by rotating a two-dimensional wedge with the opening 
angle cp around the axis perpendicular to one of the walls and intersecting the wedge's apex. 
This confining geometry might allow - under appropriate thermodynamic conditions - for a 
filling transition (23, [26]. Thermodynamic states of the fluid are parametrized by temperature 
T and chemical potential /x, and are assumed to be away from the bulk liquid-gas critical point. 

Our macroscopic analysis is based on the grand canonical functional Q([f],T, fx, h, ip) which 
is parametrized by the thermodynamic state T, fj, and geometric parameters h, (p. The allowed 
shapes of the liquid-gas interface r = f(z) are assumed to be cylindrically symmetric. The 
equilibrium shape f(z; T, fj,, h, ip) minimizes the functional f2 and leads to the thermodynamic 
grand canonical potential &(T, /j,,h,<p) = Q([f(z;T,fj,,h,ip)],T,n,h,ip). To shorten notation 
we shall refrain from displaying the dependence of analyzed functions and functionals on 
parameters T, /i, h, ip from now on. Comparison of the grand canonical potential values cor- 
responding to configuration with and without the bridge allows us to tell which configuration 
is favorable under given thermodynamic and geometric conditions, and to construct the rele- 
vant phase diagram. In what follows we shall work with the functional Afi[/] representing the 
difference between the grand potential corresponding to the state with a bridge and without 
it. For practical reasons, instead of temperature T we switch to angles 0j, i = 1,2 related to 
the temperature-dependent surface tension coefficients cjj 9 , an, and ai g , i = 1, 2 via the Young 
equation [8] 

Pig -(Til = CTlgCOsOi, £ = 1,2 , (1) 

where i stands for the £-th substrate. In what follows, the above relation will be applied also 
to thermodynamic states slightly off the bulk coexistence. 

The functional Afi[/] contains the relevant surface free energies and the appropriate bulk 
contribution for states off the bulk coexistence. Contrary to the two-dimensional version of 
the this problem [2,7], in the present case the line tension Tu g contributions related to the 
three phase contact lines where gas, liquid and the £-th substrate meet should be taken into 
account. The grand canonical functional takes the following form 
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where parameters A = oi g j ApAp with Ap = po — p > 0, Ap = p\ — p g > 0, and fu g = m g /ai g 
have dimension of length. The symbols @(z) and S(z) denote the Heaviside and Dirac function, 
respectively. 

The equilibrium interfacial shape f(z) minimizes AJ1[/] and fulfills the equation 



supplemented by two boundary conditions 
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where the coordinate z 2 is such that ffa) = 0(^2). The lhs of Eq.j3|) contains the mean 
curvature [28J and thus the surface of the bridge has constant and negative mean curvature; 
it is called a concave nodoid [29, 30J. 
Note that after introducing the contact angles 71 and 72 
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Eqs. ([4j take the form of modified Young equations [H 
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Unfortunately Eq.j3|) supplemented by the above boundary conditions, Eqs. ([!]), can not be 
solved analytically [32]. One has to resort to numerical procedures and we choose the shooting 
method [33J. For a given starting point z 2 > h the quantities f(z 2 ) and f'(z 2 ) are chosen 
in accordance with the boundary condition in Eq. (jJJ), and the function f(z) is constructed 
towards the point z = by solving Eq. ((3|). Once the point z = is reached, the boundary 
condition in Eq. j6j) is checked. If it's fulfilled then the constructed function is accepted as 
f(z) (see Fig. [2]); if it's not then the whole procedure must be repeated for a new choice 
of the starting point z 2 . The existence of the solution of Eq. (J3J) depends on the choice of 
the thermodynamic state and the values of the geometric parameters. If there is no solution 
f(z) to Eq.j3j) then the stable phase of our system corresponds to the absence of a liquid 
bridge. Fig. [2] shows interfacial shapes (broken lines) in case of identical substrates, fixed 
distance h, zero line tension coefficients, and three different choices of angle ip; for reference, 
the corresponding shapes in two-dimensional case with the same choice of parameters are 
plotted (solid lines). 
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3 Phase diagrams 



When the grand canonical potential difference is evaluated for solution f(z) of Eq.j3j) then 
depending on the sign of Afl[f] three cases are possible: (a) A$7[/] < - the bridge phase is 
favorable, (b) Afi[/] > - phase without bridge is favorable, (c) Af2[/] = - two previous 
phases coexist along a line in the (T, //)-plane which will be denoted as /i = ^afm{T). 

Before presenting the results for the currently studied three-dimensional case we recall that 
similar analysis performed for the two-dimensional case showed that both the Kelvin law for 
capillary condensation in a slit of width h and the complete filling transition in a wedge (with 
the opening angle ip) can be treated as special cases of bridge formation [27]. In particular, 
the lowest temperature at which the bridge can be formed (for given angle ip) coincides with 
the wedge filling temperature Tf(cp) (9(Tf(ip)) = 7r/2 — c^/2) which is an increasing function of 
if [25, 26] (system's geometry constrains the range of angle ip to [0,7r/2]). Upon crossing the 
coexistence curve the system exhibits discontinuous phase transition: a bridge with circular 
liquid-gas interface is formed whose width £h 1/A/z additionally depends on the opening 
angle <p and becomes infinite for parallel walls, i.e., for <p = 0. The dependence of the phase 
diagram on parameters h and A// enters via the product /iA/i. 

The phase diagrams obtained for the three-dimensional system in the case of vanishing line 
tensions, i.e., for fu g = fn g = do not differ qualitatively from those in two dimensions, see 
FigEl Both for two and three dimensions the coexistence line Hafm is a decreasing function of 
the angle ir/2 — 9. Obviously, the coexistence lines in two- and three-dimensions are identical 
in the special case of a slit corresponding to ip = 0. Similarly as in two-dimensions, the 
filling transition temperature plays a distinguished role; the AFM coexistence lines meet the 
bulk coexistence tangentially at the value of parameter ir/2 — 9 corresponding to the filling 
temperature Tf (ip). This particular value is denoted as n/2 — 9*, where 9* = 9(Tf(ip)). Fig|4] 
presents the numerically obtained plot of ir/2 — 9* as function of angle ip. It is worthwhile to 
note that the value 9* = is reached for ip = ipo < ir/2, where ipo = 1.42. Technically, ipo 
presents the largest value of angle cp for which inscribing a catenoid (corresponding to angle 
9 = 9*) into AFM-geometry is possible. 

The above results were obtained within macroscopic analysis which does not take into 
account the interaction between the liquid-gas interface and the substrates as specified e.g., 
by the interface potential uj(f) [34J, [35]. We note that the formulation of this problem on 
the mesoscopic scale in which the interface potential plays important role meets from the 
very outset some basic questions related to the structure of the effective Hamiltonian relevant 
for the present geometry in which the substrate's curvature varies along the AFM tip. One 
possible way of dealing with these problems is to start the analysis from microscopic level as 
specified, e.g. by the density functional theory [36j-[39|. 

4 The role of line tension 

In order to determine the influence of line tension on the phase diagram's structure one must 
know the line tension coefficient's dependence on the thermodynamic state of the system 
within the range of interest here. In spite of rather vast literature on this subject [40j-|44j 
there is still lack of general statements applicable outside immediate vicinity of the wetting 
points where universal behavior is observed. We recall that the line tension coefficient Tu g can 
be of either sign and its value changes significantly in the vicinity of the wetting temperature 
T w i, [45]- [47]. The behaviour of the line tension coefficient near T w i depends on the order 
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of the wetting transition and on the type of interaction between particles constituting the 
system. The universal property of this dependence is that the line tension coefficient is an 
increasing function for T /* T u ,j, and the limiting value Tu g (T w ) is positive for the first-order 
wetting (can be even infinite for specific choice of interactions), and is zero for critical wetting 

m 

In view of the lack of precise information on the dependence of line tension coefficients on 
chemical potential and temperature we shall estimate their influence on the phase diagram 
by assuming them to be constant parameters. We shall examine the system properties corre- 
sponding to different constant values and, in particular, both signs will be taken into account. 
For simplicity, in the following analysis we again consider identical substrates, i.e., 6\ = 62 = 9 
and ?u g = f 2 i g = t. 

The line tension influences the grand canonical potential not only through its presence in 
Eq.([2]) but also via boundary conditions where it affects the contact angles, Eqs (|6|7|) . The 
coexistence lines corresponding to different values of the line tension are shown on FigEl For 
positive r-values the shape of the coexistence line does not change qualitatively with respect 
to the t = case; it is still monotonic, i.e., Pafm is an increasing function of 9 and achieves 
the lowest value at the wetting temperature where 6(T = T w ) = 0. 

On the other hand, for r < the coexistence lines are not monotonic; they exhibit a local 
minimum and approach the bulk coexistence line \i = for T — > T w , see FigEl One can argue 
for this type of behaviour at small 9 values in the following way based on Eqs([6j [7J: close to 
the wetting temperature both f(z2) and /(0) must increase in order to keep cos 71, cos 72 < 1. 
Thus the size of the bridge which is proportional to the (//o — p)~ l also increases which means 
that /jlafm approaches /xq- 

The presence of line tension also influences the temperature at which the coexistence lines 
approach (tangentially) the bulk coexistence. For positive line tensions this particular tem- 
perature does not depend on the actual r-value. For negative line tensions this temperature 
decreases with | r |. This behaviour is illustrated on Figs El [3 

For two-dimensional systems [27] and for three-dimensional systems with r = there is only 
one quantity with dimension of length which can be built out of the system parameters, 
namely A = ai g / ApAp. The equation for the coexistence line (corresponding to Afi[/] = 0) 

can be rewritten as Tq (o, ^) = 0, where Tq is a dimensionless function. Thus the phase 
diagrams displayed in variables ir/2 — 9 and p — Po(T) measured in ai g /Aph do not depend 
any more on height h. If one considers system with non-zero line tension then there is one 
additional parameter with the dimension of length, namely f. This time the equation for the 
coexistence curve takes the form T T \ 9, ^, ^ =0 and the coexistence lines corresponding to 
different /i-values are not identical unless f is properly rescaled, see FigEl 

5 Sizes of liquid bridges 

The crucial aspect of the present analysis is that it deals with macroscopic objects. In Fig[2]we 
see that the shape of the bridges's surface depends on the distance from the bulk coexistence 
measured by p, — po. In this paragraph we would like to analyze the size of the liquid bridge 
in more detail, and to see how it depends on parameters describing the system. 

As the measure of the width t of the bridge we take the minimal value of f(z) in the 
range of z = to z = 22- In the case of macroscopic bridges considered in this paper their 
width must be large compared to the bulk correlation length £ which we take to be of the 
order £ pa O.lnm. The way the width I depends on height h, angle (p, and the line tension are 
depicted on Figs IMTTl For temperatures 9 < 9*((p) the width of the bridge grows to infinity 
upon approaching the bulk coexistence whereas for 9 > 9* (if) the width of the bridge achieves 
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a fixed value for \i — > /j-q(T). Obviously, the latter situation takes place only for negative line 
tension coefficients, r < 0. 

For small h- values and f « 0.1£ (for typical values of r = 10~ n J/m and ai g = 10~ 2 J/m 2 
the coefficients f = r/ai g are in the nanometer range) our macroscopic analysis predicts for 
<p close to 7r/2 the existence of bridges with width smaller than £. One has to treat these 
results with caution because such situation requires more detailed approach than the present 
macroscopic description, see [48] for similar arguments in the case of wedge filling transitions. 
We note that the minimal size of the bridge included into our analysis can not be smaller then 
the largest among the parameters h and f which have the dimension of length. This restriction 
reduces the range of parameters for which reliable phase diagrams can be constructed within 
present approach. 

6 Summary 

We analyzed macroscopically the formation of liquid bridges in AFM-like geometry in three 
dimensions with the microscope tip modeled by cone. In the case when the line tension is 
not taken into account the properties of the phase diagram displaying the coexistence of two 
phases, one with the liquid bridge present and the other without bridge, is similar to the 
phase diagram obtained for two-dimensional version of this problem. The coexistence line 
is an increasing function of the contact angle 6, and it meets the bulk coexistence line tan- 
gentially at the value of the contact angle 9*(<p) corresponding to the filling transition. We 
provide numerically obtained plot of tt/2 — 9* versus <p which points - at least for small (p 
values - at linear dependence with prefactor 1. We also note that the particular value 9* = 
corresponding to wetting is obtained not for tp = ir/2 but for smaller value ip = (po = 1.42. 
In the presence of non zero line tension values, in particular negative ones, the shape of 
the coexistence line changes substantially. For negative values of line tension coefficients the 
coexistence lines exhibit local minima and reapproach the bulk coexistence (p, = fio) for tem- 
peratures close to T w i. In addition we analyzed the size of the bridge present in the system and 
displayed its dependence on system parameters. Our investigations are limited to macroscopic 
bridges and to situations in which line tensions are treated as constants, i.e., do not depend 
on the thermodynamic state of the system. The requirement that the characteristic lengths 
of this problem are larger than the bulk correlation length reduces the range of parameters at 
which reliable conclusions can be drawn; this is depicted on Figs 9-11. 

A more detailed analysis based on the capillary Hamiltonian [25], [26] might allow one to 
estimate the role of fluctuations on the formation and stability of liquid bridges investigated 
in this paper. We also point at difficulties related to such mesoscopic formulation of the bridge 
formation problem. 
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Figure 1: A liquid-like bridge (/) surrounded by the gas phase (g) connects the planar (1) and the 
conical (2) substrates. 
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Figure 2: The shapes of bridges in two- (solid lines) and three-dimensional (broken lines) cases 
for specific choice of parameters if = ir/6, Q\ = 9 2 = n/6, fu g = f 2 i g = 0, and (a) A = lOh, (b) 
A = 50h, (c) A = 100/i. The variables r and z are displayed in h units. 



9 




Figure 3: The schematic phase diagram in variables \i — fi ( measured in aig/Aph units) and 
7r/2 — 6 displaying the coexistence lines in three- (dotted lines) and two-dimensional system (solid 
lines) of phases in which the liquid bridge is present and absent, respectively, for the case of 
identical substrates (#1 = #2), vanishing line tension coefficients (fu g = f 2 i g = 0), and for different 
values of the angle f. (a) <p = 0, (b) ip = ir/6, (c) ip = ir/4. 




Figure 4: The value of angle tt/2 — 9* at which the filling transition takes place as function of angle 
<p. The broken line represents the results of numerical analysis described in the text for the case 
r = 0; the solid line is obtained under the assumption that catenoid describes the bridge's surface. 
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Figure 5: Coexistence lines corresponding to different values of line tension: (a) fu g = r 2tg = —2h, 
(b) fn g = f 2 i g = -h, (c) fug = f 2 i g = 0, (d) fug = f 2lg = h, (e) fu g = f 2lg = 2h for specific choice 
of ip — 7r/6. The variable \i — Hq{T) is measured in aig/Aph units. 




Figure 6: The values of line tension for which the corresponding coexistence line intersects the line 
(a) A = Ai = I0 3 h (solid line), (b) A = A 2 = 10 4 /i (broken line), (c) A = A 3 = I0 5 h (dotted line). 
For A — » oo (/i = /io) the resulting plot becomes vertical at n/2 — 6* for positive r. The opening 
angle is fixed at <p = n/6. 
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Figure 7: Numerical derivatives of the coexistence curve fi = Hafm{k '/2 — 6) with respect to its 
argument denoted as fi\ FM for different line tensions (fu g = f 2 i g = f): (a) f = — 2h, (b) f = — h,(a) 
f = 0,(d) f = /i,(e) f — 2h for specific choice of <p = n/6 and 9i — 9 2 — 9. 
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Figure 8: The coexistence lines corresponding to fixed r = la.u. and for different heights: (a) 
h = la.u, (b) h = 2a. u., (c) h = 3a.u.. For r = the coexistence lines for different h are identical, 
curve (d). The opening angle is fixed and equal <p = 7r/6. The difference \i — fio(T) is measured in 
aig/Aph units. 




Figure 9: The bridge's width £ in £ units (upper graph) and the corresponding phase diagram with 
the coexistence curves (lower graph) displayed for parameter values at which macroscopic approach 
is applicable, i.e., for £ > 1 for different /i- values: (a) h = 0.1£, (b) h = 0.3£, (c) h = 0.5£, (d) 
h = 0.7£, (e) h = 0.9£. The angle cp — 7r/6 and r = 0. The difference /x — fio(T) is measured in 
aig/ApC, units. 
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Figure 10: The bridge's width I in £ units (upper graph) and the corresponding coexistence curves 
(lower graph) displayed only for parameter values at which macroscopic approach is applicable for 
different angles: (a) <p = 0.5, (b) <p = 0.7 and (c) <p = 0.9. The height h = 0.25£ and r = 0. The 
difference /i — (J,o(T) is measured in ai g / Ap£ units. 




Figure 11: The bridge's width £ in £ units (upper graph) and the corresponding coexistence curves 
(lower graph) displayed for parameter values at which macroscopic approach is applicable (b) for 
different values of line tension: (a) r = -0.2£, (b) r = -O.lf , (c) r = 0, (d) r = O.lf , (e) r = 0.2f . 
The height h = 0.1£ and y5 = 7r/6. The difference /i — Hq{T) is measured in ai g /Ap^ units. 
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